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Abstract 

The possibility of a solid state electron electric dipole moment (EDM) experiment based on 
Gadolinium Iron Garnet (GdIG) is investigated. GdIG appears to exhibit superparamagnetism 
and this effect can be used to enhance the electric-field-induced EDM signal at relatively a high 
temperatures, as compared to a simple paramagnetic system such as Gadolinium Gallium Garnet. 
The sensitivity of a GdIG based experiment might be large enough that an improvement by over 
three orders of magnitude, compared to the existing electron EDM limit, could be obtained at a 
modest temperature of 4 K, assuming that the superparamagnetic effect is not suppressed at this 
temperature, and that the spin relaxation time remains sufficiently rapid. 

PACS numbers: 
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I. INTRODUCTION 

Following the original suggestion of Shapiro we have been investigating the possi- 
bility of an improved electron EDM measurement by use of paramagnetic insulating solids 
together with modern magnetometry Q]. The basic idea is that when a paramagnetic in- 
sulating solid sample is subjected to an electric field, if the constituent atoms or ions of 
the solid have an electric dipole moment, the atom or ions spins will tend to become spin- 
polarized. Because the atoms or ions carry a magnetic moment, the sample acquires a net 
magnetization. We previously identified Gadolinium Gallium Garnet (GdGG) as possibly 
the best paramagnetic system for this study, and have done preliminary measurements in- 
dicating that the anticipated sensitivity can_be achieved, assuming that spin-glass effects do 
not enter at temperatures of order 10 mK {3J. Unfortunately, such effects appear to become 
important near 0.1 K J4], and this will limit the noise and ultimate sensitivity of a GdGG 
based experiment. We are continuing to investigate the ultimate sensitivity of a GdGG 
based experiment, however, according to the data in 4j, the magnetic susceptibility has a 
peak near 0.1 K. 

We also have performed an experiment based on the converse effect: When a ferromag- 
netic sample is magnetized, an electric polarization develops within it when the constituent 
ions or atoms carry an EDM [5] . To understand this experiment fully, the initial (zero field) 
magnetic permeability of the materials used in that work (Gadolinium Yttrium Iron Gar- 
nets) needs to be known. In a literature search for the initial permeability of these materials, 
a single reference was found (|, and an interesting effect is evident in the data presented 
there for GdlG. Shown in Fig. 1 are the values for the real part of the relative permeability 
jj! (from jj, — y! + ifi") multiplied by the measurement temperature and plotted against 
temperature. The linearity of the plot is striking, and shows a paramagnetic behavior. This 
is rare for ferrites and ferromagnetic materials; usually, the susceptibility falls rapidly with 
temperature and this fact led us to consider pure paramagnetic systems such as GdGG. An 
earlier solid state experiment that employed a nickel zinc ferrite suffered from this loss in 
susceptibility at low temperature. l| 

For solid state electron EDM experiment that employs an applied electric field, materials 
with as as large of a magnetic susceptibility as possible are desired. This is because the 
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EDM-induced magnetic field signal is (mks units) 

B = X deE*/^ a (1) 

where x is the magnetic susceptibility, d e is the electron EDM, E* is the effective electric 
field (net electric field at the lattice site of interest times the relativistic atomic or ionic 
enhancement factor) , and fi a is the magnetic moment of the paramagnetic atoms or ions that 
carry the electron EDM. We see immediately that x should be as large as possible. Note 
that the sample temperature does not explicitly enter into this formula. We are therefore 



free to use whatever means possible to increase x- As can be seen from the data in [Oj, the 
relative permeability of GdIG is fi' = x + 1 — 58 at 77 K; assuming the paramagnetic 1/T 
behavior holds for GdGG, this is already equal to the permeability of GdGG at 40 mK. 
If the trend in relative permeability shown in Fig. 1 continues to sufficiently temperature 
(4 K), a sensitive experiment, using GdIG instead of GdGG, will be possible. Scaling the 
results of Q to a relative initial permeability of 1000 (four times that of GdGG at 10 mK), 
with sample area 10 cm 2 , indicates an experimental sensitivity of 2 x 10~ 27 ecm-^sec, with 
10 kV/cm applied to the sample at a modest temperature of 6 K. Thus, in one second of 
measurement, a sensitivity at the level of the best experimental limit for the electron EDM 
8j can be achieved. 

II. PREVIOUS WORK 

There have been results from three solid state electron EDM experiments reported. The 
first of these, mentioned above, used a nickel-zinc ferrite, with \J = 11 at 4 K. Our recent 
work employed GdGG, and obtained a better value than Q], but not as large as we had 
anticipated. This is due in part to the wrong values of the electric field being used in [3j, |7|. 

The internal electric field at a site on a cubic lattice site is given by [9|] as eE / 3 where e is 
the dielectric constant, and Eq is the externally applied field strength, i.e., the field between 
the sample field plates without the dielectric material. However, the analysis presented 
in jt| is applicable only when the frequency of the applied field is much higher than the 
ion vibration frequency or principal phonon modes as there is an implicit assumption that 
the ion position does not change with application of the field. Clearly, these are optical 
frequencies, and the results must be modified for a static field. This problem was carefully 
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treated in where the change in equilibrium position due to the applied field is taken 



into account, see Sec. VI of that paper. A more refined calculation is presented in [12J, and 
when modified by a more refined relativistic enhancement factor [kJ, the energy shift of a 



Gd +3 ion, due to an electron EDM of 1.6 x 10 26 e cm (the present experimental limit [8|) 
is 

Ae = 3.8 x 1(T 22 eV (2) 

for an internal electric field of E = 10 kV/cm (externally applied field E = eE). It should 
be noted that for this calculation, e m 30 was used, whereas we found e = 19 for both single 
crystal and poly crystalline samples, and that the electron EDM energy shift is proportional 
to a factor (see [111], Eq. (44)), 

x = fr-ljg = fe-ljg, 1 

327renGd e 327renGd 

and in the limit of large e, the sensitivity is independent of e. We can calculate the apparent 
electric field at a Gd +3 ion site, by taking the ionic enhancement factor as 2.2 from [13| |. 

3.8 x 10- 22 eV = 2 2 x 1£ x 1Q _ 27 ecm ^ent ^ Kpparent = 3 6 w/cm (4) 
oil £/0 

when Eq = 10 kV/cm, so the field is reduced by a factor of three at the ion site compared 
to the applied field. Interestingly, this is the field expected within a spherical cavity in a 
dielectric material, for large values of dielectric constant. Clearly, the value E = eE /3 used 
in [?J is probably not correct, although the correction in that case is relatively small because 
e = 2.2 for the nickel-zinc ferrite used in that work. 

Furthermore, when using a different value for the dielectric constant for different GdGG 
samples, the factor X discussed above must be modified. As mentioned already, this factor 
makes the experimental sensitivity independent of e for e >> 1. This means that the limit 
obtained in relation to the efforts [3j are a factor of 1.5 too small, taking e w 20 compared 



to e = 30 in 



111 ], without modifying X. Also, a density correction has to be applied also 



that varies between samples. 



III. SUPERPARAMAGNETISM: GADOLINIUM IRON GARNET 



Superparamagnetism exists in a ferromagnetic system when the size of the domains be- 
comes so small that the interaction energy between domains, or the energy required for 
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domain wall movement, is below the thermal fluctuation energy, h{I '. A magnetized ferro- 
magnet is actually in a metastable state, and the relaxation time to the ground state for a 
domain-domain interactions can be described by an Ahrrenius-type equation, 



T = r exp(K u V/k b T) (5) 

where K u is the uniaxial anistropy constant, and V is the volume of a particle or domain, 
[lij] In this context, it is useful to think of the spin of a domain as undergoing a Brownian 
motion around the easy axis, with K u being the energy to change alignment from one easy 
axis to another. Alternatively, we can apply the same equation to the movement of a domain 
wall (e.g., increase in domain volume), with K U V replaced by an appropriate "activation 
energy" associated with an increase in domain volume, in which case 

t = r exp(E d /k b T) (6) 

where Ed is the energy barrier for domain wall movement. It is generally accepted that at low 
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frequencies and low fields, the relaxation time is determined by domain wall movement. 

If the time (e.g., the inverse of the frequency of an oscillating magnetic field) of an ex- 
perimental measurement of a magnetic property is longer than r, a ferromagnetic material 
can behave as a superparamagnetic system, with the properties that the susceptibility in- 
creases as 1/T and that there is very little magnetic hysteresis. On the other hand, for 
faster times, hysteresis becomes important, and the system appears as a usual ferromagnet. 
Normally, V is large, and r becomes infinite for all practical purposes. Nevertheless, there 
are well-known instances where V is small enough for r to be small, the most famous of 
which is the superparamagnetic limit on the density of a hard drive of a few hundred gigabits 
per square inch. A smaller bit size, as determined by the thickness of the magnetic layer 
(to restrict the transverse field leakage distance) is unstable against thermal fluctuations. 
Increasing this density is of course an area of active research. Although not named as such, 



superparamagnetism was first investigated by Oliver Heaviside in 1887 [15] . He showed that 
finely powdered iron dust, "Worked up with wax into solid cores (1 wax to 5 or 6 iron by 
bulk)," showed no hysteresis effects, or, in his words, the effect of the iron in such a state 
is, "Exactly equivalent to self-inductance." 

For the EDM problem, superparamagnetism is useful, as we have discussed already. As 
shown in Fig. 1, we see that the real part of the relative initial permeability of GdIG can 



be described by 

, /mN 5830 K 
fi'(T) = — 17.8 (7) 

where the constant part probably represents the finite accuracy of the data, but does suggest 
that /i' = 1 at 310 K, which is well below the Curie temperature of GdIG, 590 K, but is close 
to the "compensation temperature" where the three sublattice magnetization cancel [9J. A 
question remains as to the range of validity of this relation; it is probably not overly hopeful 
to assume that it is unmodified to temperatures approaching 4 K. It is interesting to note 
that in jf|, for GdYIG materials, the relative permeability falls with temperature. It is of 
course not possible to tell if this is due to the domain size, due to a different value of K u , 
or due to domain movement. In [6J], the "dispersion frequency" (peak in fi" as a function of 
frequency) increases for GdYIG mixtures at the temperature is lowered; this is the opposite 
of what is expected for superparamagnetism. There are no dispersion frequencies listed 
for GdIG except for 77 K where the first dispersion frequency is given as 27 MHz. It is 
impossible to know how the scaling varies with temperature based on this single value; in 
addition, it is likely that only a small subset of the sample spins contribute to the observed 
effect. Different domains can have different Ej,, so as the temperature is lowered, there will 
likely be a gradual decrease in r. 

As the correlation time decreases, the noise due to spin fluctuations appear in a narrower 
frequency range, and therefore have a larger contribution near zero frequency, where an 
EDM experiment would be performed [13]. In a fundamental sense, if the relative permeabil- 
ity continues to increase as 1/T, at sufficiently low temperature, the correlation time will 
increase as exp(EdfkbT) as will the zero frequency spin noise, and will eventually outpace 
the increase in fx' = 1 + y\ the real part of the susceptibility. 

Using the results in [19j, we can determine the effects of the relaxation time in a quan- 
titative fashion. The density of Gd +3 in GdIG is pad = 1-24 x 10 22 /cc, and the magnetic 
moment comes mostly from the Gd which has fi a = 7 fib where fib is the Bohr magneton. A 
signal to noise of unity indicates the level with which an EDM can be discriminated above 



background noise. The results in 19J can be recast to give a fundamental EDM sensitivity, 
in terms of the sample volume V, 

, Liofx a 2y/2p Gd Ti . . 

"e = ~ /= (8) 



XeE*\/V 

so to achieve d e = 1.6 x 10~ 27 ecm in one second of averaging, for E* = lOkV/cm, V — 10 



cc, and x — 1000 requires that T\ < 10 10 sec, or a relaxation frequency of about 1 GHz. 
On the other hand, we can ask what level of sensitivity can be achieved at 77 K, with \ = 57 
and 7\ = 1/(2tt x 27MHz) = 5 x 10~ 9 sec and E* = 10 kV/cm: d e = 2 x 10~ 25 eony/sec, so 
after ten days of averaging, a limit of 2 x 10 -28 ecm can be achieved. Therefore a substantial 
improvement in the experimental limit appears as imminently feasible if even cooling below 
50 K or so is not possible. 

Of course, GdIG can be permanently magnetized when a large enough magnetic field is 
applied. The arguments above assume that the applied fields are very small. For GdIG, 
the super- paramagnetism is a non-linear effect, in that it disappears above some critical 
field, of value to be determined. An earlier study of the susceptibility of GdIG has a large 
applied field so the superparamagnetic effect is not evident. 16] More recent studies of 
the magnetic properties of GdIG show that the susceptibility remains relatively large at 5 K. 
However, the slope, hence inferred susceptibility, at zero field in Fig. 9 of 17J is determined 
by the sample geometry; for a spherical sample in a solenoid, the effect of the "demagnetizing 
field" limits the apparent susceptibility to x = M/ Bo = 3 in the limit of large x ( or A 4 ) ( see 
Eq. (5.115) of 18]), which is close to inferred value of x f° r an of the curves in Fig. 9. 

The maximum field before a sample might be magnetized significantly and permanently 
can be estimated by analogy with usual ferromagnetic materials, which have paramagnetic 
properties at low fields while in the demagnetized state. The sample begins to acquire a 
significant magnetization in an external field when the applied field is larger than (neglecting 
the geometry-dependent demagnetizing field) 

Bo = ^ (9) 

X 

where M s is the saturation magnetization. If we take x = 1000, M s = p a pGdi we find that 
B = 10 Gauss. 



IV. POSSIBLE EXPERIMENT 



As shown in Fig. 2, operation of a solid state EDM experiment at temperature near 4 K 
will allow the use of a "warm" -bore system as heat loads are not so much of an issue. The 
heat shields and dewar will have to be made of materials with low electrical conductivity 
(plastics, G-10), and their thermal conductivity could be improved by incorporating thin 
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wires which would not contribute as much Johnson noise as continuous metal sheets. The 
entire system could be easily contained in a superconducting shield. 

The magnetic field due to an EDM of d e = 1.6 x 10~ 27 ecm, with E* = 10 kV/cm, 
X = 1000 is 

B = x^l w 0.05 fr. (io) 

Ha 

By placing a superconducting pickup loop with inductance L p around the sample, the flux 
$ = BA = L p I can be converted to a current. This current can be measured in two ways, 
either by placing another loop around magnetometer cell, or by sending the current to a 
SQUID magnetometer. 

For the first option, the present best magnetometer noise has been achieve in a spin- 
exchange- relaxation- free (SERF) potassium (K) magnetometer, at the level of 0.1 fT/\/Hz 
[20! ] . If the sample is relatively large compared to the magnetometer, it will be possible 
to have a smaller loop around the magnetometer and step up the field at the magnetome- 
ter. The optimal transfer of signal occurs when the magnetometer coil inductance matches 
the pickup coil inductance. An order of magnitude increase in the field magnitude at the 
magnetometer is not beyond comprehension. 

A simpler approach is to use a SQUID magnetometer. If we treat the SQUID as a 
current measuring device, its sensitivity is determined by the fractional flux quantum noise, 
typically l/i^o/v^Hz for commercial grade SQUIDs, and an order of magnitude better for 
research units. The input current is converted to a SQUID flux through a coupling coil's 
mutual inductance, which is typically 0.2 /xA/$o [2l| . where $0 is the superconductor flux 
quantum. Taking a sample area of 10 cm 2 , using B — 0.1 fT, and Lp = 5 /zH (which 
is about a factor of two larger than the SQUID input inductance), a current I p = 20 fA 
is generated, implying a SQUID signal of 0.1 /i$o/"\/Hz, a factor of 10 below the intrinsic 
SQUID sensitivity . It is unlikely that use of a transformer to better match L p to the input 
inductance would help as most any pickup loop will be within an order of magnitude of the 



SQUID input inductance. 
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We see immediately that the magnetometer method offers at least a factor of 10 better 
signal to noise that the SQUID method. The magnetometer method can be further im- 
proved by use of a flux concentrator loop, as mentioned above. The possibility of such an 
improvement when using a conventional magnetometer was already mentioned in [2j. On 
the other hand, the simplicity of the SQUID system is quite attractive, and the sensitivity 
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of the magnetometer configuration is possibly better than the fluctuation noise of the sam- 
ple. If this latter point proves correct, the added complexity of an atomic magnetometer is 
unwarranted. 

An even simpler experimental concept is shown in Fig. 3. Here, two disks of GdIG 
sandwich a magnetometer, and although the loss in signal is at least a factor of five for 1 cm 
thick disks, separated by 3 cm, the sensitivity is comparable to the SQUID method outlined 
above. 

V. REQUIRED STUDIES 

There are a number of issues that need to be addressed to evaluate the efficacy of GdIG 
for a solid state electron EDM experiment. The following is a listing, not necessarily in 
order of importance or complete, of points to be addressed: 

1. Determine the effective electric field E* for GdIG which includes both the atomic en- 
hancement and the dielectric shielding. There is no reason to suspect that there should be 
surprises compared to previous calculations of GdGG. 

2. Calculate the effective electric field susceptibility, noting that the electric field couples 
significantly to only the Gd +3 spins. Most of the magnetization of GdIG comes from Gd +3 
at temperatures below 77 K. 

3. Verify the 1/T scaling of the permeability to low temperature, and measure fi" as a 
function of frequency and temperature. The latter should give the most sensitive measure 
of the domain relaxation time, and will provide a means to determine the minimum oper- 
ating temperature. These are the most important parameters in regard to the design of an 
experiment. It could be that only a small fraction of the domains contribute to fi" so the 
scaling in [l9| must be modified; the value of T% is no longer simply determined. Also the 
issue of domain rotation vs. growth needs to be studied. If domain rotation is the principal 
contribution to the relaxation time, then the results in [3] again need to be modified. 

4. Establish the maximum magnetic field that can be applied before the sample begins to 
acquire a permanent magnetization. This is related to the low-frequency hysteresis, which 
is important regarding permanent magnetization due to displacement, currents, microdis- 
charges, and leakage currents. 

5. Determine the effects of sample preparation on the superparamagnetic properties. As the 
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crystalline structure of the sample becomes coarser, it might be expected that the super- 
paramagnetic effect will be suppressed, due to the V factor in the exponent of the Arrhenius 
time, or an increase in the energy required for domain wall movement. 
6. Determine the degree to which residual sample magnetization can be eliminated. The 
limiting residual magnetism will place limitations on the types of magnetometry that can 
be used. Also the stability to the residual magnetization is important. Given that the Curie 
temperature of GdIG is 590K, it should be possible to anneal a sample while it is placed in 
a magnetic shield and heated within an oven. In addition, the compensastion temmperature 
of GdIG is near room temperature, so careful handling and avoidance of exposure to large 
magnetic fields would be necessary to ensure that the residual magnetization remains low. 
A crucial related point is addressed in (4) above. 

VI. CONCLUSION 

GdIG exhibits the behavior normally associated with a superparamagnetic system, so 
it is probably fair to refer to GdIG as having this property. Unfortunately, the initial 
permeabilities of the rare earth garnets have not been extensively studied, and a significant 
amount of effort will be required to assess the suitability of GdIG for a solid state electron 
EDM experiment. Nevertheless, GdIG has potential to provide the basis of an experiment 
that can achieve a truly significant improvement in the electron EDM limit. If it is possible 
to achieve // « 1000 for GdIG by cooling to 4 K, after ten days of averaging, an EDM 
limit at the level of 10 -29 ecm can be achieved using a commercial SQUID magnetometer. 
Furthermore, GdIG will be useful as a cryogenic ferrite and will certainly find applications in 
cryogenic electronic systems. In particular, if the permeability remains large at 4 K, GdIG 
will be very useful in the construction of low noise amplifiers, including DC parametric 
amplifiers. 

As a final note, there are possible applications to other fundamental measurements. By 
searching for a GdIG magnetization correlated with velocity or direction in absolute space, 
Lorentz violating interactions can be sought, with perhaps to six orders of magnitude im- 
provement over existing limits. 
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VIII. APPENDIX: DETUNING OF THE G RADIOMETER CONDITION DUE 
TO SAMPLE PERMEABILITY 



In our work[3], we employed a superconducting pickup loop, operated in a concentric 
loop planar gradiometer configuration. Unfortunately, the permeability of the sample alters 
the pickup loop inductances, and spoils the gradiometer tuning, and reduces the experiment 
sensitivity slightly because the loop inductance is increased. In the case of a single loop or 
radius a around an infinitely long cylinder of permeability fi and radius b, the increase in 



inductance is 



23] 



/■oo 

SL = 2^ a2 / (jfea)]2tZjfe (11) 

Jo 

where /x is the permeability of free space, Ki(ka) is the first order modified Bessel function, 
and 

- l)kbh{kb)h{kb) ^ I (kb) 
{ ' (ji- l)kbK (kb)h(kb) + 1 ~ K (kb) 1 ' 

where Jo and K are the zero-order modified Bessel functions, and the approximation holds 

when fj, » 1. The effect of finite length is small, provided that length i > 2b, which is the 

case for the experiment. In the case a = b for a loop directly on a cylinder, 5L « 6/io« ~ 

0.1 /iH, for a = 1.2 cm, a 20% effect. When a = V2b, 5L = 0.02 /iH, a much smaller effect. 

The loss in sensitivity is around 10%, leaving the principal problem due to this effect as 

the detuning of the gradiometer condition, which is temperature dependent, and of order 

10%. Usually, gradiometers are fabricated to provide 1% common mode rejection. We see 

a substantial limitation to the degree of common mode rejection possible when there is a 

permeable sample in the central loop. 



[1] F.L. Shapiro, Sov. Phys. Uspekhi 95, 145 (1968). 
[2] S.K. Lamoreaux, Phys. Rev. A (2002). 

[3] C.-Y. Liu and S.K. Lamoreaux, Mod. Phys. Lett A 19, 1235 (2004). 

11 



[4] P. Schiffer et al., Phys. Rev. Lett. 74, 2379 (1995). 

[5] B.J. Heindenreich et al., Phys. Rev. Lett. 95, 253004 (2005). 

[6] G.E. McDuffie, Jr., J.R. Cunningham, Jr., and E.E. Anderson, Jour. App. Phys., Supplement 

to Vol. 31, 47S (1960). 
[7] B.V. Vasil'ev and E.V. Kolycheva, Sov. Phys. JETP 47, 243 (1978). 

[8] B.C. Regan, E.D. Commins, C.J. Schmidt, and D. DeMille, Phys. Rev. Lett. 88, 071805 
(2002). 

[9] Charles Kittel, Introduction to Solid State Physics, 4th Ed. (Wiley, New York, 1971). 
[10] Richard M. Bozorth, Ferromagnetism (Van Nostrand, Princeton, N.J., 1951). See, e.g., Chap. 
18., where it is pointed out that the usual mechanism for magnetization is by domain move- 
ment. 

[11] S. Kuenzi, O.P Sushkov, V.A. Dzuba, and J.M. Cadogan, Phys. Rev. A 66, 032111 (2002). 
[12] T.N. Muckhamedjanov, V.A. Dzuba, and O.P. Sushkov, Phys. Rev. A 68, 042103 (2003). 
[13] V.A. Dzuba, O.P. Sushkov, W.R. Johnson, and U.I. Safranova, Phys. Rev. A 66, 032105 
(2002). 

[14] L. Neel, Ann. Geophys. C.N.R.S. 5, 99 (1949). 

[15] O. Heaviside, The Electrician, Sept. 25, 1885, p. 375. (Reprinted in Electrical Papers by Oliver 

Heaviside, vol. II, p. 159 (Macmillan and Co., London, 1892).) 
[16] W.P Wolf and R.M. Bozorth, Phys. Rev. 124, 449 (1961). 
[17] T. Yamagishi et al., Philosophical Magazine, 85, 1819 (2005). 
[18] J.D. Jackson, Classical Electrodynamics, 2nd Ed. (Wiley, New York, 1975). 
[19] D. Budker et al., Phys. Rev. A 73, 022107 (2006). 
[20] J.C. Allred et al, Phys. Rev. Lett. 89, 130801 (2002). 
[21] Quantum Design, Technical Specification for DC SQUID. 
[22] Quantum Design, SQUID Application Note 1052-202. 

[23] W.R. Smythe, Static and Dynamic Electricity (McGraw-Hill, New York, 1950); Chap. 8, prob. 
15. 



12 




50 



100 



150 200 
T [K] 



250 



300 



FIG. 1: The values of GdIG relative initial permeability taken from 
T. The linear form shows that [/,' oc 1/T+ constant. 
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FIG. 2: Schematic of a cryostat with a warm bore that incorporates a SERF magnetometer. Two 
disks of GdIG have opposite sign high voltage applied to electrodes on the disk faces, with a 
ground plane iplaced between the disks. A SQUID magnetometer could be incorporated instead, 
in which case a second warm bore is not required. If it turns out that the magnetic properties of 
GdIG remain favorable to 4 K (or below), the warm bores can be dispensed with if a SQUID is 
employed. However, as pointed out in the text, better sensitivity can be achieved with an atomic 
magnetometer, in which case a single warm bore would be required, with the GdIG held at 4 K. 
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FIG. 3: Schematic of a cryostat with a warm bore that incorporates a SERF magnetometer between 
two disks of GdIG that have opposite sign high voltage applied and a ground plane is placed between 
the disks. For disk separation larger than the thickness of the disks, there is a substantial loss in 
magnetic field signal at the magnetometer location, compared to the internal field of the GdIG 
samples. 
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